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Abstract 

A formula is derived for stiffness of a polymer chain in terms of the distribution function 
of end-to-end vectors. This relationship is applied to calculate the stiffness of Gaussian chains 
(neutral and carrying electric charges at the ends), chains modeled as self-avoiding random 
walks, as well as semi-flexible (worm-like and Dirac) chains. The effects of persistence length 
and Bjerrum's length on the chain stiffness are analyzed numerically. An explicit expression is 
developed for the radial distribution function of a chain with the maximum stiffness. 
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1 Introduction 

This paper is concerned with stiffness of polymer chains and explicit formulas for its determina- 
tion. As this issue lies on the border between statistical physics of macromolecules and mechanical 
engineering, it has not attracted substantial attention in the past. The situation has changed dra- 
matically in the past decade due to the development of novel experimental techniques for the direct 
measurements of force-extension relations on individual chains in bio-polymers. A correct measure 
of stiffness is important for our understanding of the differences between the responses of flexible 
and semi- flexible chains, as well as for the assessment of the effects of segment interactions (for 
example, excluded-volume interactions for neutral macromolecules and electrostatic interactions 
for polyelectrolyte chains) on their mechanical properties (see ^ and the references therein). 

In the statistical physics of polymers, the stiffness of a macromolecule is conventionally associ- 
ated with the so-called "effective spring constant:" the ratio of a force applied to a free end of a 
chain (its other end is assumed to be fixed) to the displacement of the free end along the force direc- 
tion An advantage of this method is that the stiffness is determined in terms of the response of 
an individual chain, which is convenient from the experimental standpoint. A shortcoming of this 
approach is that it may lead to conclusions that contradict the physical intuition. As an example, 
we refer to the fact that this "macro" -stiffness decreases with an increase in the "micro" -stiffness 
characterized by the bending rigidity of a semi-flexible chain. A detailed discussion of this issue is 
provided in Section 2. 

Mechanics of polymers focuses on the response of ensembles of macromolecules and treats an 
elastic modulus of an ensemble as a natural measure of its stiffness. A (shear or Young's) modulus 
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is expressed in terms of the strain energy density of a polymer, and it reflects the cumulative 
contribution of the strain energies of chains and the energies of their interaction. When the inter- 
chain interactions are negligible (dilute polymer solutions) or their energy may be adequately 
described in a simple way (the influence of mutual interactions between chains in rubbery polymers 
and polymer melts is traditionally accounted for in terms of the incompressibility condition [3J, an 
elastic modulus of an ensemble divided by the number of chains per unit volume provides a natural 
measure of stiffness for an individual macromolecule. 

We apply the latter approach and derive an explicit expression for the stiffness of a chain 
"embedded" into some ensemble subjected to affine deformations. The objective of this study is 
to demonstrate that the "mechanical" definition of stiffness of a macromolecule is free from the 
above shortcoming (an increase in the bending rigidity of a semi-flexible chain induces the growth 
of its stiffness, as it may be expected), on the one hand, and that it results in the same formula 
for the stiffness of a Gaussian chain that the conventional definition of stiffness does, on the other. 
To develop an analytical formula for the stiffness, we calculate the strain energy of an ensemble 
of chains under an arbitrary affine deformation, find the Young's modulus at uniaxial tension, and 
associate the stiffness of an individual chain with its contribution into the modulus. This procedure 
allows an explicit expression to be derived for the stiffness of a polymer chain with an arbitrary 
end-to-end distribution function. The formula is illustrated by several examples, where the stiffness 
is expressed in terms of parameters of the distribution function, and the dependence of stiffness on 
these parameters is analyzed numerically. 

The exposition is organized as follows. In Section 2, the conventional stiffness of a chain fi is 
calculated, and it is found that /i is inversely proportional to the mean square end-to-end distance 
of a chain. The strain energy of an ensemble of polymer chains with an arbitrary end-to-end 
distribution function is determined in Section 3. The formula for the strain energy is simplified for 
uniaxial tension with small strains in Section 4. An explicit expression for the stiffness of a chain 
is derived in Section 5. Section 6 focuses on several examples, where the stiffness is determined 
analytically, and the effect of parameters of the distribution function on this quantity is studies 
numerically. Some concluding remarks are formulated in Section 7. 

2 Conventional stiffness of a chain 

We begin with the conventional in statistical physics definition of stiffness for a polymer chain. Our 
aim is to develop an explicit expression for this parameter and to demonstrate that it is inversely 
proportional to where b is the mean square end-to-end distance of the chain. 

A chain is treated as a curve with some length L in a three-dimensional space. An arbitrary 
configuration of the chain is described by the function r(s), where r stands for the radius vector 
and s G [0, L]. For definiteness, we assume the end s = to be fixed at the origin, r(0) = 0, and 
the end s = L to be free. An "internal structure" of the chain is characterized by a segment length 
bo and a number of segments N ^ 1, which are connected with L by the formula L = boN. 

A chain is entirely determined by its Hamiltonian H{r{s)). Given H, the distribution of end- 
to-end vectors Q is described by the propagator (Green's function) 



where is Boltmann's constant, T is the absolute temperature, and the path integral is calculated 
over all curves r(s) that start at the origin and finish at the point Q, 




(1) 



r(0) = 



r(L) = Q. 



(2) 
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The path integral in Eq. (1) is determined unambiguously by the normalization condition 



G{Q)dQ 



For a chain loaded by a force F at the free end, the Green function reads 



Gf{Q) 



r(L)=Q 
r(0)=0 



exp 



1 



Hir)-F-Q) V[ris)] 



(3) 



(4) 



where the last term in the exponent stands for the work of external force. As this term is indepen- 
dent of the curve r(s), Eqs. (1) and (4) imply that 



Gi.(Q) =G(Q)exp 



F Q 



V 



Substitution of expression (5) into the partition function Zp = J Gi?(Q)(iQ results in 



Zf = / G(Q) exp 



F Q 



dQ. 



Differentiation of Eq. (6) with respect to F implies that 

kuT dZp 



(Q) 



Zp dF ' 



(5) 



(6) 



(7) 



where (Q) = / QGF{Q)dQ J GF{Q)dQ stands for the average end-to-end vector. Formula (7) 
means that the vectors (Q) and F are collinear, while their moduli are connected by 



(Q) 



k-BT dZp 
Zf dF ' 



At small forces F, Eq. (8) can be linearized, 



(8) 



(9) 



and the chain stiffness is described by the effective spring constant ^ in Eq. (9). 

To reveal a disadvantage of this definition, we consider the Marco-Siggia interpolation formula 
1^ for the force-stretch relation of a semi-flexible (worm-like) chain, 



k^Tyl 



Ir 



1 u 
4 + 1. 



-4V L' 

where Ip stands for the persistence length. It follows from this equality and Eq. (9) that 

SksT 



1^ 



2/pL 



(10) 



Bearing in mind that the persistence length Ip of a semi-flexible chain is connected with its bending 
rigidity k by the formula Ip = KL/{k-QT)^ we find from Eq. (10) that the stiffness ^. is inversely 
proportional to the bending rigidity n. This conclusion appears to be counter-intuitive, because 
one expects that the larger the rigidity of a chain is at the micro-level (described by the parameter 
k), the higher its stiffness ^ is at the macro-scale jS]. 
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To show that this shortcoming is independent of the concrete form of the force-stretch relation, 
we derive an exphcit expression for the pre-factor /x in Eq. (9) confining ourselves to chains with 
isotropic Green functions G = G.^{Q), where Q = |Q|. Introducing a spherical coordinate frame 
{Q, (j), 6} whose z vector is directed along F, we find from Eq. (6) that 

roo r2-rr r-rr / ms ^ \ 

Zf= G^{Q)Q^dQ d(t) exp( 7 ^ jsingdg. 

Jo Jo Jo ^ KB-L ' 

Calculation of the integrals over ^ and Q results in 

Zp = 47r^ J^^ G,(Q) sinh(0;)QdQ. 

Substitution of this expression into Eq. (8) implies that 

+ Io°"G4Q)^oM-^)Q'dQ 
F /-G,(Q)sinh(g)QdQ ■ 

Setting Q = bx, where the mean square end-to-end distance b reads 



J^G^{Q)Q^dQ 

and introducing the notation 



f = ^. u = % (12) 



we arrive at the force-stretch relation 

„ "'^a) 1 ,,,, 

with 

poo poo 

Wi{f)= G^{bx)smh{fx)xdx, W2{f)= G^bx) cosh{fx)x'^dx. 
Jo Jo 

At small dimcnsionlcss forces /, the hyperbolic functions are expanded into the Taylor series in /, 

which implies that 

roc ^3 /-oo 

Wlif) = f G^{bx)x'^dx + ^- G4bx)xUx + . . . , 
Jo Jo 

roc ^2 i-oo 

W2if) = / G4bx)x^dx + — / G^{bx)x^dx + . . . , 
Jo 2 Jo 

where the dots stand for terms of higher order of smallness. Substituting these expressions into 
Eq. (13), neglecting terms beyond the first order of smallness, returning to the initial notation, 
and using Eq. (11), we find that 

/ = 3n (/ « 1). 
This equality together with Eq. (12) yields 

F = ^{Q) m«b). (14) 

Equations (9) and (14) imply that 

H = 3A;B^r^ (15) 
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which means that the conventional stiffness /x of a chain with an arbitrary end-to-end distribution 
function is a merely geometrical parameter that is inversely proportional to the square of the 
average end-to-end distance b. According to Eq. (15), the stiffness of a Gaussian chain with the 
mean square end-to-end distance be reads 

liG = 3A:bT65'. (16) 

Obviously, Eq. (16) coincides with Eq. (10) with the persistence length Ip = ^bQ. Setting ^' = 
h/HGj we find from Eqs. (15) and (16) that 



3 Strain energy of a chain in an ensemble 



Our aim now to derive a formula for the strain energy of a chain in an ensemble of macromolecules 
whose deformation at the micro-level coincides with macro-deformation (the affinity hypothesis). To 
simplify the analysis, we assume the deformation to be incompressible and adopt the conventional 
hypothesis that inter-chain interactions may be accounted for by the incompressibility condition. 

Denote by Q the end-to-end vector of a chain in the initial (reference) state and by q the 
end-to-end vector in the actual (deformed) state at an arbitrary instant t > 0. The distribution 
functions of end-to-end vectors in the initial and deformed states read po(Q) ^-nd p{t,q), respec- 
tively. Transformation of the reference state of the chain into its deformed state is described by 



q = F(t)-Q, 



(18) 



where F is the deformation gradient for macro-deformation, and the dot denotes inner product. 
The function F(t) obeys the differential equation 



F(0) = I, 



(19) 



where L(t) is the velocity gradient, and I is the unit tensor. For an incompressible macro- 
deformation, the Smoluchowski equation for the function p{t, q) reads jSj 



dp dp 



L-q, p(0,q) =po(q)- 



dt c?q 
The solution of Eq. (20) is given by 

p{t,ci)=po{F'\t)-(i), 
and it satisfies the normalization condition 



J p(*,q)c?q = 1- 



(20) 



(21) 



(22) 



The distribution functions Po{Q,) and p{t, q) are expressed in terms of appropriate configurational 
free energies C/o(Q) a-nd U{t,q) by the Boltzmann relations 



PoiQ) = exp 



pit, q) = exp 



U{t,q) 



(23) 
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There are two ways to determine the strain energy of a chain W. According to the first, we calculate 
the increment AU of the configurational free energy caused by transition from the reference state 
to the actual state, 

AU{t, Q) = U{t, Q) - f/o(Q) = -kBT[lnp{t, Q) - Inpo(Q) 
and average it with the help of the distribution function in the reference state, 

Wi{t) = -/cbT J [lnp{t, Q) - lnpo(Q)]po(Q)dQ- (24) 

Following the other approach, the increment of the configurational free energy is calculated with 
respect to the actual state, 

AU{t, q) = C/o(q) - U{t, q) = -A;Br[lnpo(q) - lnp{t, q 
and it is averaged by using the distribution function in the deformed state, 

W2{t) = -ksT J [inpo(q) - lnp(t, q)]p(t, q)dq. (25) 
Substitution of Eq. (21) into Eq. (24) results in 

Wi{t) = -kBT J [lnpo(F-i(t) • Q) - lnpo{Q)\po{Q)dQ. (26) 
Combining Eqs. (21) and (25), we find that 

W2{t) = kBT J [\npo{F-\t) ■ q) - lnpo{ci)\po{F-\t) ■ q)dq. 

Introducing the variable Q by Eq. (18) and bearing in mind that dq = dQ for an incompressible 
deformation, we arrive at the formula 



W2it) = kBT / Inpo(Q) - lnpo(F(t) • Q) Po{Q)dQ. 



(27) 



It seems natural to define the strain energy of a chain W as the weighted sum of the strain energies 
Wi and W2 calculated by using different ways of averaging of the configurational free energy, 

W ={1- a)Wi + aW2, 

where a G [0, 1] is a material parameter. Substitution of Eqs. (26) and (27) into this formula 
implies that 



W{t) = kBT J a(lnpo(Q)-lnpo(F(0-Q)) + (l-a)(lnpo(Q)-lnpo(F-^(i)-Q)) 
For an isotropic distribution of end-to-end vectors in the reference state, we set 

Po(Q) = P{Q^), 



PoiQ)dQ. (28) 



(29) 



where = Q ■ Q, and P{r) is a given function of a scalar argument r. Combining Eqs. (28) and 
(29) and taking into account that 

(f • Q) • (f • q) = Q • C • Q, (f-1 • q) • (f-1 • q) = Q • • Q, 
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where the left and right Cauchy-Green deformation tensors read 



B = F F^, C = F^ F, 

and T stands for transpose, we find that 

W = kBT J a(lnP(Q2) _ lnP(Q . C • Q)) + (1 - a)(lnP{Q^ - lnP(Q • B"^ • Q)) 
It follows from this equality and the identity 

JlnP{Q ■ B-^ • Q)P(Q2)c^q = J\nP{Q ■ ■ Q)P{Q^)dQ 

that the strain energy per chain reads 



(30) 



P{Q^)dCl. 



W = kBTj a(lnP(Q2)_inP(Q.C-Q)) + (l-a)(lnP(Q2)_inP(Q.C-^-Q)) 



P{Q^)dQ. (31) 



Neglecting the energy of inter-chain interaction, we calculate the strain energy per unit volume of 
an ensemble of chains W as the sum of the strain energies of individual chains. 



W = kBTM J a(lnP(Q2)_inP(Q.C-Q)) + (l-a)(lnP(Q2)-lnP(Q-C-^-Q)) 

where M is the number of chains per unit volume. 



P{Q^)dQ, (32) 



4 Elastic modulus of a chain 

Our aim now is to apply Eq. (32) in order to determine the elastic modulus of an ensemble of 
chains under uniaxial tension 



xi = kXi, X2 = k 2X2, X3 = k 2X3, 



(33) 



where k is an elongation ratio, {Xi} are Cartesian coordinates in the reference state, and {xj} are 
Cartesian coordinates in the deformed state {i = 1,2,3). It follows from Eqs. (18), (30) and (33) 
that 

C = fc^eiei + /c~"^(e2e2 + 6363), C"^ = fc^^eiei + /e(e2e2 + 6363), (34) 

where ej are base vectors of the Cartesian frame {Xj}. Substituting expressions (34) into Eq. (32) 
and introducing spherical coordinates {QjCpjO}, whose z axis coincides with ei, we find that 



W 



k^TM 



00 r2-K 

P{Q^)Q^dQ / , 

JO 



a ( In P{Q')- In P(Q2 (^2 ^^^z q ^ ^-i ^:^^z g 



1,-1 „;„2 , 



+(1 - a) (in P(Q2) _ InP(Q2(A;-2 cos^ 9 + /csin^ 9))) 



smOde. 



Performing integration over (f) and introducing the new variable x = cos 9, we obtain 

W roo _ _ rl 



kBTM = ^T^l"p{Q^)Q^dQj^^a(lnP{Q^)-lnP{Q^{ex^ + k-\l-x^)))) 
+ (1 - a)(lnP(g2) _ lnP(Q2(yfc-2^2 ^ _ ^2 



dx. 



(35) 
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Equation (35) provides an exact formula for the strain energy density of an ensemble of chains 
under uniaxial tension. Simple algebra (see Appendix) implies that at small strains, when 

k = l + e (e<l), (36) 

this equation reads 

A7re\Ki + K2), (37) 



knTM 

where terms beyond the second order of smallncss are disregarded, and the coefficients Ki and K2 
are given by (the prime denotes the differentiation) 



= - / P'{Q')QUQ, K2 = - 

Jo 5 Jo 



P(Q2 



P"{Q 



2^ 



Q'^dQ. (38) 



It is worth noting that the dependence of W on the parameter a disappears at small strains. It 
follows from Eq. (29) that for an isotropic Green function G, 

P{Q') = G.iQ), P'{Q') = ^, ^"(^') = i^[^*^^)-^]- (39) 
Equations (38) and (39) imply that 

Ki = -\ GUQ)Q'dQ, K2 = ^ ^^^QUQ - ^ j^iG:{Q)Q^ - G:{Q)Q^)dQ. 
Bearing in mind that 

poo roo 

/ G'i{Q)Q*dQ = -4 / G',{Q)Q''dQ, 
Jo Jo 

we find that 

Substitution of these expressions into Eq. (37) results in 



W = ^TTkBTMe' r ^'^Q. (40) 

Jo Lr*{Q) 

At uniaxial tension of an incompressible medium, the strain energy per unit volume W reads 

W = Ee^, (41) 

where E denotes Young's modulus. It is worth noting that Eq. (41) does not contain the coefficient 
^ on the right-hand side that conventionally arises at small uniaxial deformations. It follows from 
Eqs. (40) and (41) that 

^=?.fe™f (42) 
For an ensemble of Gaussian chains with the end-to-end distribution function 



«-(« = (24)'-K-i|)- («) 

Eq. (42) reads 

18 / 3 \ - /""^ / 30^ \ 

Ea = ^.k,TM{^yl exp(-^)QVQ. 

Setting z = QV^/bc and calculating the integral, we find that 

Eg = ^kBTM. (44) 
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5 Stiffness of a polymer chain 



Our aim now is to determine stiffness of a polymer chain based on Eq. (42) , which coincides (to some 
extent) with Eq. (9). The remark in the parentheses refers to the fact that definition (9), widely 
used in the statistical physics of macromolecules, differs from that in mechanical engineering where 
the stiffness is measured as the ratio of an appropriate force to a strain (not to a displacement). 
To avoid this discrepancy, we replace {Q) in Eq. (9) by the ratio {Q)/b, which, together with Eq. 
(16), implies that the conventional stiffness of a Gaussian chain is given by 



(45) 



To obtain a counterpart of Eq. (45) grounded on Eq. (42), it seems natural to introduce a 
hypothetical ensemble of closely packed chains (the latter means that the number of chains per 
unit volume M equals v~^, where v = |vr6'^ is the average volume occupied by a chain), and to 
define the stiffness of a chain S as the Young's modulus of this ensemble. If follows from Eq. (44) 
that for a Gaussian chain, 

which differs from Eq. (45) by an insignificant pre- factor (of order of unity) only. In the general 
case, this definition implies that 



S 



3kBT 



Q^dQ. 



1063 Jo G4Q) 
Excluding the coefficient b with the help of Eq. (11), we arrive at 



S 



3kBT 
10 



G,{Q) 



Q^dQ 



!^G,{Q)Q^dQ 



I^G,iQ)Q*dQ 



(47) 



(48) 



Formula (48) allows the stiffness of a chain with an arbitrary end-to-end distribution function 
G*((5) to be calculated. Introducing the relative stiffness $ as the ratio of the chain stiffness to 
that for a Gaussian chain, <I> = S/Sq, we find from Eqs. (17), (46) and (47) that 



Ait 3 

^2 

15 



G,{Q) 



Q^dQ. 



(49) 



6 Examples 

Our purpose now is to calculate the ratio <I> for several distribution functions G^{Q) and to show 
that Eq. (49) leads to physically plausible dependencies of on material parameters. 



6.1 A flexible chain with excluded- volume interactions 

The end-to-end distribution function of a fiexible chain modeled as a self-avoiding random walk 
may be approximated by the stretched exponential function [HIIZI, 



(50) 



G*{Q) = gexp 

Here 6 and / are positive constants, and the pre-factor g is found from the normalization condition 



1 

47r 



exp 



■) )Q^dQ 



2TrPT{ 



25 > 
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where r(x) = exp(—z)z^ ^dz is the Euler gamma function. Substitution of expression (50) into 
Eq. (48) imphes that 



5 



O/cbT 
407r/3 



(3 + 26) 



r( 



25' 



(51) 



When 6 = 1, we set / = ^cy § in accord with Eqs. (43) and (50). In this case, Eq. (51) is reduced 
to Eq. (46) for the stiffness of a Gaussian chain. It fohows from Eqs. (46) and (51) that 



2(5 + 3 



3r(i) 
2r(i) 



(52) 



Setting 5 ^ in Eq. (52), we find that ^ 
brackets. Bearing in mind that 



due to the decay of the expression in the square 



hm 

5— too 



nil 



we conclude that 



hm — 

5— ►oo 



5 
3' 

5 
2' 



which means that $ increases hnearly with 6 at sufficiently large values of the exponent. To 
analyze the effect of 6 on the relative stiffness of a chain, we calculate ^ for S S [1,2], i.e. in the 
interval where typical values of this parameter are located. The results of numerical simulation are 
presented in Figure 1, which shows that the stiffness monotonically increases (practically linearly) 
with 6. As d may be treated as a measure of strength of excluded-volume interactions, this implies 
that repulsive segment interactions cause the growth of chain stiffness. It is worth mentioning that 
our conclusion contradicts Eq. (17): according to the latter formula, an increase in 6 causes the 
growth of the average end-to-end distance b, which, in turn, induces a decrease in the conventional 
dimensionless stiffness ^. 



6.2 A Gaussian chain with electric charges at the ends 

Our aim now is to assess the effect of an electrostatic field on the stiffness of a Gaussian chain with 
two equal charges e fixed at its ends. It is assumed that the mean square end-to-end distance of the 
chain be is smaller than the Debye screening length / d , which implies that the Coulomb interaction 
between charges is not screened (Hj. The Hamiltonian of a non-charged Gaussian chain reads 



2bo Jo ^ds 

To account for the energy of interaction between the charges, we replace Hq by the Hamiltonian 

where e is the dielectric constant of an ion- free dilute solvent in which the chain is immersed. 
The last term on the right-hand side of Eq. (53) describes the energy of electrostatic repulsion of 
charges. Substituting Eq. (53) into Eq. (1) and bearing in mind that the last term is independent 
of the curve r(s), we find that 



G(Q)=.exp(-||)exp(-|). (54) 
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The first exponent in Eq. (54) stands for the Green function for a neutral Gaussian chain, 
e^/{ekBT) is the Bjerrum length, and the constant g is determined by Eq. (3), 



with = iBV^/bc- Formula (54) correctly predicts that the probability to find the end s = L of 
the chain in the close vicinity of the origin (the position of the other end s = 0) strongly decreases 
due to repulsion of charges. It follows from Eqs. (46), (48) and (54) that 



^ roc / ^2 5 



5/o^cxp(-f ^ ^Jz^dz L/o°°exp(-f - ^)z^dzl 



/-exp(-4-|)z2d,il 



(55) 



The function ^{(,) is plotted (in the double logarithmic coordinates with log = log^g) in Figure 
2. According to this figure, $ decreases with ^. At sufficiently large ^, the curve ^(0 ^^-y be 
approximated by the dependence 

log$ = $0 - ^ilogC, (56) 

where the coefficients $o and $i are determined by the least-squares method. Our results of 
numerical analysis demonstrate that $i ^, which implies the scaling law 

$oc(^)^ (/b»M- (57) 

The conclusion that the stiffness of a charged Gaussian chain is smaller than that of an appropriate 
neutral chain may be explained by the fact that electrostatic repulsion of chain ends increases the 
end-to-end distance in the reference state, which means that the number of available configurations, 
and, as a consequence, the chain entropy are reduced. It is worth noting, however, that the same 
explanation is inapplicable to flexible chains with excluded-volume interactions, because Eq. (52) 
demonstrates the growth of stiffness due to segment interactions. 

6.3 A flexible chain with the maximum stiffness 

It is of interest to determine the distribution function of end-to-end vectors for a chain whose 
stiffness is maximal. The analysis is confined to smooth radial distribution functions G*((5), which 
are positive at any Q G [0, oo), tend to zero rather rapidly as Q — ^ oo, and remain bounded together 
with their derivatives in the vicinity of Q = 0. We suppose also that the second moment of the 
distribution function (Q^) = b'^ and the most probable end-to-end distance Ic arc fixed. The latter 
quantity is determined from the condition of maximum of the function G^{Q)Q^, 

Ic = argmaxG*((5)(5^. 
Q 

To demonstrate that these conditions uniquely determine the radial distribution function of a chain 
with the maximum stiffness, we present Eq. (47) in the form 



where the functional reads 



iv'iQ)?^..^ . ......... (59^ 



11 



The presence of the last term in Eq. (59) allows arbitrary (non-normalized) Green functions y{Q) 
to be considered. Denote by A the maximum of the functional J^{y) on the set of smooth positive 
functions y{Q). A function yo{Q) that maximizes J^{y) satisfies the equality 



f 

Jo 



yoiQ) 



dQ = 0. 



(60) 



The first term in Eq. (60) is transformed by integration by parts 



f 

Jo 



yoiQ) 



-dQ 



y'oiQ)Q' 
yoiQ) 



dyoiQ) = yo{Q)Q' 







yoiQ) 



dQ. 



Our assumptions regarding the properties of the function yo (Q) imply that the out-of-integral term 
vanishes. Combining this equality with Eq. (60), we obtain 



r\( y'oiQ)Q^ \' , 
Jo 



yoiQ) 



yoiQ)dQ = 0. 



This equation is fulfilled provided that the function yoiQ) obeys the differential equation 



^ yoiQ) ' 



Integration of Eq. (61) implies that 



y'oiQ) 



0. 



B 



A 



yoiQ) Q^ 3Q' 

where B is an arbitrary constant. The general solution of Eq. (62) reads 



yQiQ) = AQ 3exp(-^), 

/here A is another constant. The parameters A, B and A arc found from the equalities 

47r / yoiQ)Q''dQ = 1, Att yoiQ)Q^dQ = hi, yo(Q)Q' ^ = 0, 

JO Jo '- -iy=<c 



(61) 



(62) 



(63) 



(64) 



which describe the normalization condition for the distribution function, the definition of the mean 
square end-to-end distance, and the definition of the most probable end-to-end distance, respec- 
tively. Substitution of Eq. (63) into the last formula in Eq. (64) results in 



B 



A- 6 3 
3 



(65) 



Dividing the second equality in Eq. (64) by the first and using Eqs. (63) and (65), we find that 



/o°° exp( 



(A-6)fg 



)Q-^dQ 



/o-exp(-(^)Q-^dQ 



Setting z = Q/lc and calculating the integrals, we arrive at the formula 

A-6^§r(A^) 



r(¥) 



(66) 
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For a Gaussian chain with be = be and Ic = ^Gy |; set A = 9Ai and find that Ai obeys the 
transcendental equation 

A 2N| r(Ai-f) 3 



It fohows from Eqs. (46) and (58) that 



' (68) 



'S'max 27Al 



Solving Eq. (67) for Ai numerically and using Eq. (68), we obtain Ai = 3.1809 and Sq = 
0.0582 5max- The result is rather surprising. It means that Gaussian chains are not so flexible: 
the stiffness of a Gaussian chain is about 6% of the maximal stiffness of a chain with the same 
geometrical parameters. 

The difference between the shapes of the radial distribution functions for a Gaussian chain and 
for a chain with the maximum stiffness is seen in Figure 3, where the results of numerical simulation 
are presented for bo = 1.0. The shape of a Gaussian chain is described by Eq. (43), whereas the 
shape of a chain with the maximum stiffness is determined by Eq. (63), where B is given by Eq. 
(65), A is calculated from Eq. (67), and A is found from Eq. (3). Figure 3 demonstrates that the 
distribution function of end-to-end vectors for a chain with the maximum stiffness vanishes in the 
vicinity of the point Q = 0, has a pronounced maximum at Q = Ic, and slowly decreases with Q 
when Q > Ic- 

It is worth noting the importance of the assumption regarding the smoothness of the distribution 
function G^{Q) and its positiveness in [0,oo). If these constraints are violated, a distribution 
function of end-to-end vectors may be constructed for a chain with an infinite stiffness. An example 
is given by 

G*(Q) = Afl-f?)'! W<0, G*(Q)=0 (Q>/). (69) 



Formula (69) and similar expressions for the radial distribution function naturally arise for polymer 
chains treated as random walks governed by equations of anomalous diffusion ^QlEj. It should be 
emphasized, however, that Eq. (47) is inapplicable to chains, whose distribution functions have 
finite supports. This is explained by the fact that Eq. (47) is grounded on the hypothesis regarding 
an affine deformation of a network, while the probability p(t, q) is not defined by Eq. (20), when 
the vector q = F • Q is located in the domain where the distribution function po(Q) vanishes. 

6.4 A semi-flexible chain 

Semi-flexible chains provide another class of macromolecules, for which the applicability of Eqs. 

(47) and (48) may be questioned. This is explained by the fact that the integral in Eq. (32) for the 
strain energy density diverges. On the other hand, appropriate integrals in Eq. (48) converge, and 
it is tempting to employ this formula in order to calculate the stiffness. Our aim now is to assess 
the effect of persistence length 1^ on the stiffness of worm-like chains and to demonstrate that Eq. 

(48) results in a physically plausible behavior of the function S{lp). 

Assuming that G^{Q) = A{Q/L), where L stands for the chain length, and setting z = Q/L, 
we find from Eq. (48) that 



?,k^T {k{z)f 4 



S = / — . , , z dz 

10 Jo K{z) 



S^k{z)z^dz 



!^k{z)z^dz 



(70) 
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We begin with the analysis of semi-flexible chains with the distribution function 



A(z) = (z > 1). (71) 

Here C = Ip/L stands for the dimensionless persistence length, H2{z) = Az^ — 2, and the pre-factor 
g is determined by Eq. (3). We substitute expression (71) into Eq. (70), calculate the integrals 
numerically (by the Simpson method with the step Az = 0.001) taking into account 500 terms in 
the series, and calculate S. The ratio S = S/{k^T) is plotted versus C in Figure 4 at relatively 
small values of C (in the linear scale) and in Figure 5 at arbitrary Q (in the double logarithmic 
scale) . 

The results of numerical simulation based on Eq. (71) are compared with those found by using 
the radial distribution functions proposed in [T2l [T3]. 



A(^) = — ^exp-— — ^— (z<l), (72) 



(1 - z2) 



2 



9 



8C(l-z2) 



A(^) = : 17#^ — exp- ^^^" (z<l). (73) 



(l-z2)2(2-z2)3 



3 



4C(l-z2) 



Equations (72) and (73) presume the function A(z) to vanish at z > 1, in accord with the second 
equality in Eq. (71). The ratios S calculated from Eqs. (70), (72) and (73) are plotted versus the 
dimensionless persistence length C, in Figure 5. The results of numerical analysis are approximated 
by the function 

log 5 = ^0 + Si log C, (74) 

where the coefficients So and Si are determined by the least-squares technique. Figure 5 shows 
that Eq. (74) correctly fits the dependence S{C,) at relatively large values of C- The coefficient 
Si in Eq. (74) is close to two for all models under investigation, which implies that the stiffness 
of a worm-like chain grows with persistence length /p being proportional to 1^ (in other words, 
being proportional to the square of the bending stiffness k), as it is expected 0. However, some 
discrepancies between the scaling prediction 

5 oc (75) 

with a = 2 and the numerical results should be mentioned. The exponent a practically equals 
2 for distribution function (72), and it is close to | for functions (71) and (73) (no changes in a 
are observed with the growth of the number of terms in the series and a decrease in the step Az 
employed in numerical integration). The deviations of the exponent for the distribution functions 
(71) and (73) from a = 2 may be explained by the fact that inadequate approximations of the Green 
function were chosen in the derivation of these relations. The method suggested in is based 
on the hypothesis that any curve r(s) in Eq. (1) is close enough to the straight line connecting 
the points r(0) = and r(L) = Q. This assumption is correct for end-to-end vectors Q with 
Q = L, but it may cause large discrepancies at Q < L driven by the constraint on the curve length. 
Figure 5 reveals that the method of softening this constraint proposed in ^2] (the local restriction 
on the length of the tangent vector is replaced by the global one) leads to a substantially better 
approximation, whereas the approach developed in ^Hl does not improve substantially the quality 
of approximation. 
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Based on another way of thinking, the so-called Dirac chains were introduced in with the 
radial distribution function 

K{z) = -jJ=h[ ^'^^ ) (z<l), A(z) = (z>l). (76) 

Here Q is the dimensionless persistence length, 

is the modified Bessel function, and the pre- factor g is determined by Eq. (3). To assess the 
stiffness of a chain with the distribution function (76), we substitute this expression into Eq. (70), 
calculate the stiffness 5, and plot the ratio S versus C, in Figure 6 together with its approximation 
by Eq. (74). This figure shows that at relatively small persistence lengths, /p < 0.3-L, the stiffness 
S increases with l^, but the scaling exponent a « 1.5 is lower than that for Eqs. (71) to (73). 
With a further increase in the persistence length Zp, the stiffness diminishes, which indicates that 
the model becomes inappropriate at relatively large values of /p . 



7 Concluding remarks 

A "mechanical" definition of stiffness is introduced for a polymer chain, and an explicit formula 

is derived to express the stiffness in terms of the radial distribution function. According to our 

approach, a polymer chain is "embedded" into an affine ensemble, and the stiffness is defined as 

an elastic modulus per chain of the ensemble. 

This relation is applied to calculate the stiffness of (i) a Gaussian chain, (ii) a flexible chain 

modeled as a self-avoiding random walk, and (iii) a Gaussian chain carrying electric charges at its 

ends. The influence of the Bjerrum length 1-q on the stiffness S" of a charged Gaussian chain has 

_ 1 

been studied numerically. It is revealed that S decreases with 1-q being proportional to ^ . 

An analytical formula is derived for the distribution function of a chain with the maximum 
stiffness. It is found that the stiffness of a Gaussian chain is about 6% of that for a chain with the 
maximum stiffness and the same geometrical parameters. 

The effect of persistence length /p on the stiffness of semi-flexible chains has been evaluated 
numerically. It is demonstrated that the stiffness of worm-like chains increases with 1^ following 
the pattern E (x 1^ with a « 2. For a Dirac chain, stiffness increases with persistence length when 
Ip is less than 30% of the contour length L, and the scaling exponent a is close to |. 



Appendix 

Neglecting terms beyond the second order of smallness with respect to e, we find from Eq. (36) 
that 

A:V + k-\l -x^) = l- e(l - 3^^) + e^, 
k-'^x'^ + k{l -x'^) = l + e(l - 3x2) _^ ^^2^2^ 

The expansion of the function lnP{Q^(l + a)) in a Taylor series with respect to a in the vicinity 
of the point q = reads 

\nP{Q [l + a)) = lnP{Q ) + aQ p^giy + 2 ^ p2(Q2) + •••' (^"2) 
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where the prime stands for the derivative of the function P. It fohows from Eqs. (A-1) and (A-2) 
that 



lnP(Q2) 



In P(Q2) _ 



\uP{Q\ex^ + k-\l - X-))) = eQ\l - Sx^)^ - ^'Q'^ 



lnP(Q2(fc-2^2 ^ ^(-L _ ^2))) ^ _^g2(i _ 3^2^^"(Q') 0,2^2^2 



- Se^Q^x' 



pm 



Integrating these expressions with respect to x and taking into account that 

[\l - 3x'^)dx = 0, [\x'^dx = l, [^{l-3x'^fdx = ^, 
Jo Jo Jo 



we find that 



/ [lnP(Q2) _ lnP(Q2(fc2^2 ^ ^-i^-^ _ ^2^^^ 
Jo L 

= / \lnP{Q'^)-\nP{Q^{k-^x'^ + k{l-x'^))) 
Jo L 



dx 
dx 



p'm ^2 , HP'mr-p"iQ')p{Q') ^4 

P{Q'')^ 5 P2(Q2) ^ 

Equations (37) and (38) follow from Eqs. (35) and (A-3). 



(A-3) 
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List of figures 



Figure 1: The dimensionless stiffness $ of a flexible chain with excluded- volume interactions 

versus the dimensionless parameter 6. 

Figure 2: The dimensionless stiffness $ of a Gaussian chain with electric changes at the ends 
versus the dimensionless Bjerrum length ^. Circles: results of numerical simulation. Solid line: 
their approximation by Eq. (56) with $o = 0.031 and $i = 0.335. 

Figure 3: The radial distribution functions G*(z) with z = Q/bc for a Gaussian chain (solid line) 
and for a flexible chain with the maximal stiffness and the same parameters 6c and Ic (circles). 

Figure 4: The ratio S versus the dimensionless persistence length ( for approximation (71) of the 
Green function. 

Figure 5: The ratio S versus the dimensionless persistence length ( for three approximations of 
the distribution function. Symbols: results of numerical simulation. Unfilled circles: Eq. (71). 
Filed circles: Eq. (72). Diamonds: Eq. (73). Solid lines: their approximations by Eq. (74). Curve 
1: So = 1.27, Si = 1.83. Curve 2: Sq = 0.91, Si = 1.99. Curve 3: Sq = 0.26, Si = 1.79. 

Figure 6: The ratio S versus the dimensionless persistence length for a Dirac chain. Circles: 
results of numerical simulation. Solid lines: their approximations by Eq. (74). Curve 1: Sq = 
-0.04, Si = 1.47. Curve 2: Sq = -2.87, Si = -3.97. 
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